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NOMENCLATURE 

a  Local  sonic  velocity 

p  Local  pressure 

u  Local  flow  velocity  in  x-direction 

V  Local  flow  velocity  in  y-direction 
P  Riemann  variable 

X  Spatial  coordinate 

y  Spatial  coordinate 

t  Time  coordinate 

C  'Physical'  characteristic  for  1-D  unsteady  flow 

r  'State'  characteristic  for  1-D  unsteady  flow 

3 

p  Density,  kg/m 

Y  Ratio  of  specific  heats  for  air 

Subscripts : 

±  Representing  'families'  of  characteristics 

L  Left 

R  Right 
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I.   INTRODUCTION 

In  looking  for  30%  improvement  in  range  for  cruise  mis- 
siles, DARPA  has  recently  considered  several  alternate  con- 
cepts; the  "wave  rotor"  is  the  key  component  in  one  of 
these.   "Wave  rotors",  "wave  engines",  "wave  pressure  ex- 
changers", "wave  equalizers" — use  wave  propagation  through 
the  fluid  trapped  in  a  rotor  passage  to  transfer  energy  from 
one  fluid  to  another,  or  between  the  fluid  and  the  rotor 
shaft.   The  flow  into  and  out  of  the  rotor  is  steady  in  time 
at  ports  which  are  at  precise  locations  over  parts  of  the 
rotor  annulus.   The  rotor  passages  may  be  axial  or  helical, 
and  behave  similarly  to  shock  tubes. 

Such  "unsteady  energy  exchangers"  are  more  complicated 
than  traditional  steady-flow  gas  turbomachinery  and  hence 
have  received  less  attention.   They  are  being  reexamined 
now  because  they  offer  a  unique  advantage — they  can  be  self- 
cooling,  since  both  hot  high  pressure  gas  and  cold  low  pres- 
sure gas  use  the  same  passages  for  alternate  periods  of  time 
in  the  cycle.   Also,  the  inherently  poorer  performance  of 
steady  flow  turbomachinery  in  sm.all  sizes  has  tended  to 
limit  the  performance  of  small  gas  turbine  engines. 

In  August  1981,  Cdr .  Don  Finch  of  DARPA  asked  the  Turbo- 
propulsion  Laboratory  (TPL)  for  an  opinion  on  the  potential 
of  the  wave  engine  concept  and  positive  recommendations  were 


given  (Ref .  1) .   The  recommendations  were  followed  by  a  pro- 
posal, at  DARPA's  request,  to  begin  a  research  program  to 
develop  the  analytical  tools  to  assess  and  to  design  wave 
rotor  devices  (Ref.  2).   The  Office  of  Naval  Research  (Dr. 
A.  D.  Wood)  agreed  to  monitor  a  joint  DARP A/CNR  program 
through  a  modification  of  the  work  statement  of  the  ONR- 
sponsored  axial  compressor  flowfields  project  (Ref.  3). 

The  proposal  was  to  form  a  team  consisting  of  Dr.  Atul 
Mathur  as  full-time  investigator,  with  part-time  participa- 
tion by  Dr.  S.  Eidelman  (NRC  Associate  at  NFS) ,  Professor 
J.  R.  Erwin  (NAVAIR  Visiting  Research  Professor  in  Aeronau- 
tics), Dr.  R.  P.  Shreeve  (Director,  TPL)  and  Professor  M.  F. 
Platzer  (Department  of  Aeronautics) .   The  aim  was  to  lay  out 
an  experimental  and  analytical  program  which  would  lead  to 
a  proper  understanding  of  the  potential  and  problems  of  wave 
rotors.   The  availability  at  TPL  of  an  expensive  wave  rotor 
apparatus  formerly  used  by  Klapproth  (Ref.  4)  at  General 
Electric  and  contributed  by  Professor  Erwin,  and  its  com- 
patibility with  existing  power  supplies  and  high  response 
instrumentation,  made  it  possible  to  propose  a  basic  experi- 
mental program  at  very  low  cost.   Funding  received  June  5, 
198  2  from  ONR  allowed  a  study  program  but  not  an  experimental 
program  to  be  initiated.*   The  present  report  documents  the 


*Partial  support  by  Naval  Air  Systems  Command  and  by  the 
NRC  Associateship  program  for  different  members  of  the  team 
is  acknowledaed. 


findings  nade  by  the  TPL  team  in  the  study  of  wave  rotors  to 
31st  October  1982. 

The  work  has  concentrated  on  developing  an  understanding 
of  wave  rotors  and  of  the  methods  needed  and  constraints  in- 
volved in  their  design  and  analysis.   No  effort  has  been 
given  to  cycle  analysis  or  to  further  justifying  the  study 
of  the  technology.   Rather,  the  promise  of  the  potential 
engine  applications  is  accepted  as  proven.   This  was  the 
major  conclusion  drawn  from  a  workshop  meeting  on  wave  rotor 
technology  held  at  Mathematical  Sciences  Northwest  (MSNW)  in 
Seattle,  Washington,  6-10  August  1982  (Ref .  5) ,  which  was 
attended  by  three  members  of  the  TPL  team.   Furthermore,  a 
concentrated  study  of  the  potential  application  to  cruise 
missile  engines  is  the  purpose  of  the  current  program  at 
MSNW,  under  DARPA  sponsorship. 

Therefore,  in  the  present  report,  only  a  brief  review 
of  the  evidence  for  useful  engine  applications  is  given  in 
the  following  Section  II.   In  Section  III,  account  is  given 
of  the  tools  and  understanding  of  wave  rotor  technology  that 
have  been  developed  so  far.   A  staged  experiment  is  proposed 
which  would  provide  a  vehicle  to  evaluate  the  utility  of 
the  analytical  and  computational  models.   Conclusions  and 
recommendations  are  summarized  in  Section  IV. 


II.   POTENTIAL  FOR  ENGINE  APPLICATIONS 

In  principle,  wave  rotors  offer  the  potential  for  in- 
creasing gas  turbine  cycle  temperatixres  without  the  usual 
penalty  associated  with  having  to  cool  the  turbine  structure. 
It  is  pertinent  to  ask,  to  what  extent  has  the  practicability 
of  using  wave  rotors  in  engines  been  demonstrated?   What  has 
been  the  total  practical  experience  to  date?  We  turn  first 
to  published  reviews. 

Two  useful  reviews  of  wave  rotor  devices  have  been  writ- 
ten (Ref.  6,  7).   In  1965,  Azoury  (Ref.  6)  reviewed  the 
basic  ideas  involved  in  the  wave  rotor  concept  including  the 
fundamental  (simplified)  wave  processes  involved  in  the 
construction  of  a  useful  cycle,  wave  diagrams  (including  the 
port  opening  and  closing  processes) ,  high  and  low  pressure 
scavenge  processes,  and  a  discussion  of  the  basic  4-port 
arrangement  and  effect  of  adding  additional  ports.   Azoury 
also  reviewed  "some  industrial  applications",  including  the 
gas  generator  (wave-engine),  the  supercharger,  air  compres- 
sor, refrigerator  and  the  "equalizer/divider".   References 
were  given  for  several  experimental  programs  known  to  Azoury. 
The  references  omitted  the  work  of  Kantrowitz  (Ref.  8) ,  Von 
Chain  (Ref.  9)  ,  Pearson  (Ref.  10)  ,  and  Klapproth  (Ref.  4)  . 

The  more  recent  review  by  Rose  (Ref.  7)  acknowledged 
Azoury 's  earlier  paper,  again  reviewed  basic  ideas  and 
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updated  the  review  of  potential  applications.   Thermal  power 
cycle  applications  and  applications  in  chemical  processing 
were  dealt  with  more  extensively.   Small  gas  turbines  were 
mentioned  only  briefly,  and  no  new  examples  of  experience 
were  quoted.   Kantrowitz's  basic  work.  (Ref.  8)  was  noted, 
but  those  of  Von  Chain  (Ref.  9) ,  Pearson  (Ref.  10)  and 
Klapproth  (Ref.  4)  were  not. 

The  history  of  wave  rotor  engine  experience  can  be 
summarized  as  follows: 

Brown-Boveri  (1941-43).   Based  on  ideas  in  Seippel's 
patents,  a  direct-flow  wave  rotor  was  incorporated  in  a 
locomotive  gas  turbine  (Ref.  6).   Because  the  cycle  did  not 
give  the  promised  overall  efficiency  it  was  revised  to  use 
a  heat-exchanger  and  the  wave  rotor  was  removed.   Brown- 
Boveri  has  since  developed  and  is  marketing  the  COMPREX 
supercharger,  but  has  not  pursued  the  direct-flow  wave  rotor 
further.   It  should  be  noted  in  passing  that  the  elements 
which  are  critical  in  the  design  of  superchargers  for  diesel 
engines,  which  result  in  the  selection  of  the  reversed-f low 
arrangement  for  example,  are  somewhat  different  from  those 
which  would  be  present  in  the  design  of  a  self-cooling  gas 
turbine  wave  rotor. 

Von  Chain  (April  194  6)  .   In  Ref.  9  tests  of  a  wave  en- 
gine designed  for  1.4  kg/sec  at  full  RPM  are  reported  to 
have  oeen  made  in  which 

.  .  .  the  maximum  pressure  ratio  attained  amounted  to  2.4 
at  a  combustion-chamber  temperature  of  2200 °F.   The  engine 
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was  able  to  drive  itself;  the  effective  output  was  prac- 
tically zero;  the  scavenging  was  not  sufficient  to  remove 
all  the  hot-air  remnants  from  the  cells.  .  .  .  tests  were 
limited  to  3000  RPM  (or  45%  design  RPM)  .  .  .  because  of 
a  deficient  welding  of  the  cell  walls  to  the  rotor. 

The  historical  background  given  in  Supplement  A  of  Ref. 
9  is  reproduced  here  in  Appendix  A. 

Kantrowitz  (1949) .   A  comparison  of  experimental  mea- 
surements and  prediction  of  the  behavior  of  the  flow  in  a 
wave  rotor  were  reported  in  Ref.  8.   Kantrowitz  also  ini- 
tiated an  effort  at  NACA  Langley  which  was  continued  by  Paul 
W.  Huber  (Ref.  11) .   An  engine  was  built  and  developed 
through  several  stages  before  competing  priorities  terminated 
the  project.   No  publication  v/as  issued  and  the  apparatus 
was  sent  on  loan  to  G.E.   Nett  power  was  produced  at  "very 
low  efficiency"  (Ref.  11).   The  arrangement  of  the  engine, 
which  used  a  rotor  with  curved  gas  passages,  and  was  fed 
from  a  compressed  air  source,  is  shown  in  Appendix  B. 

Pearson  (195?-196Q) .   Reported  at  the  August  1982  Seattle 
workshop  and  in  Ref.  10,  Pearson  built  a  wave  engine  which 
output  45  horsepower  at  a  cycle  efficiency  of  9%.   This  pre- 
viously unpublished  work,  carried  out  as  a  corporately 
funded  development  project  by  Rus ton-Hornby  in  Great  Britain, 
demonstrated  that  the  prototype  wave  engine  had  superior 
part-speed  performance  to  a  contemporary  standard  gas  tur- 
bine engine  of  comparable  power  level  and,  further,  demon- 
strated that  wave  rotors  could  be  made  to  work  at  metal 
temperatures  well  below  the  high  inlet  gas  temperature. 
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The  project  showed  that  the  level  of  knowledge  necessary  to 
simply  make  a  wave  engine  work  had  been  easily  surpassed. 

Klapproth  C -1962)  .   A  research  and  development  pro- 
gram was  carried  out  at  G.S.  to  develop  a  commercial  wave 
engine.   The  most  recent  published  account  is  given  in  Ref. 
4.   It  is  noted  that  the  major  problems  were  considered  to 
be  mechanical  rather  than  gas  dynamic;  rotors  failed  when 
they  were  unshrouded  and  problems  were  encountered  with  end 
seals  due  to  circumferential  distortion.   They  did  achieve 
self-sustaining  operation  with  nett  power  output  before 
management  terminated  the  project.   The  termination  was  the 
result  of  a  policy  to  concentrate  only  on  large  engines. 
The  future  seen  for  the  wave  engine  was  in  smaller  sizes 
(Ref.  12). 

General  Pov/er  Corporation  [GPC]  (1969-19B21 .   Founded 
in  1969  explicitly  to  develop  wave  engines,  this  is  current- 
ly the  only  known  on-going  commercial  effort.   Following  the 
patent  filed  in  1971  (Ref.  13),  the  project  has  proceeded 
with  great  difficulty.   The  engine  design  is  conceptually 
sophisticated  and  employs  very  short  rotor  passages.   As  for 
any  wave  rotor,  success  depends  critically  on  whether  the 
flow  processes  have  been  correctly  calculated.   Until  opera- 
tion close  to  design  speed  can  be  sustained  (for  other  than 
structural  reasons)  the  adequacy  of  the  gas  dynamic  descrip- 
tion can  not  be  tested. 
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In  reviewing  what  has  been  achieved  toward  demonstrating 
the  practical  use  of  wave  rotors  in  engines^  the  work  of 
Pearson  stands  out  clearly  from  the  rest.   His  is  the  only 
work  which  resulted  in  a  properly  documented  measurement  of 
output  power  and  cycle  efficiency.   (It  is  entirely  regret- 
table that  the  documentation  has  only  now  become  available.) 
Most  important  perhaps  is  the  fact  that  the  machine  worked 
when  it  was  first  switched  on,  demonstrating  that  the  analy- 
tical approach  used  in  the  design  was  accurate  enough  to  get 
a  prototype  machine  working.   Like  Von  Chain  CRef.  9),  Huber 
(Ref.  11),  Klapproth  (Ref.  12)  and  MSNW  CRef.  14),  Pearson 
found  sealing  to  be  the  major  problem  (and  in  his  case,  he 
felt  the  only  problem)  standing  in  the  way  of  achieving  the 
design  performance  goal.   This  should  be  a  matter  of  some 
concern  to  GPC ,  since  the  shortness  of  their  rotor  will 
serve  to  accent  the  leakage  problem. 

The  only  other  performance  figure  for  a  wave  rotor  de- 
vice which  can  be  quoted  is  the  "energy  exchange  efficiency" 
of  70%  achieved  by  MSNW  (Ref.  14)  in  a  "tailored"  cold  ex- 
periment.  It  is  noted  that  the  geometry  in  this  case,  very 
long  rotor  passages  compared  to  the  area  of  cross  section 
and  mean  radius  to  the  passage  center  also  very  large  com- 
pared to  passage  size,  were  chosen  to  satisfy  the  assump- 
tions underlying  the  analytical  approach.   None  of  the 
engine  programs  have  used  parameters  in  a  similar  range. 
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While  the  potential  of  wave  engines  has  been  demon- 
strated definitively  by  Pearson,  the  total  documented  prac- 
tical experience  is  collectively  small.   Each  attempt  has 
been  an  experience  in  a  different  range  of  the  available 
parameters.   The  greatest  effort  by  far,  if  the  Comprex  work 
of  Brown-Boveri  is  included,  has  been  with  axial  rotor  pas- 
sages.  This  is  certainly  the  geometry  in  which  the  one 
dimensional  wave  propagation  approach  should  apply  most  ac- 
curately.  Even  constant  area,  staggered  (helical)  passages 
present  curvature  to  the  relative  flow,  and  a  calculational 
approach  which  works  for  large  rotors  with  nearly  axial  pas- 
sages may  not  be  good  enough  for  rotors  of  small  diameter 
with  highly  staggered  passages.   The  results  obtained  by 
Huber  using  cambered  passages  would  have  been  of  interest 
for  this  reason.   Pearson's  passages  were  staggered  but  at 
a  modest  angle  (-30°  to  the  axis).   The  final  Klapproth 
rotor  had  "straight"  passages  which  were  staggered  at  60° 
to  the  axis,  which  would  provide  a  more  severe  test  of  the 
purely  one  dimensional  approach. 

Thus  the  state-of-the-art  is  that  the  wave  engine  con- 
cept works,  that  a  prototype  engine  similar  to  Pearson's 
in  design  could  be  made  to  work  now  with  reasonable  risk, 
but  that  little  to  no  knowledge  is  available  to  say  which 
range  of  any  parameter  is  likely  to  give  the  best  results. 
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III.   WAVE  ROTOR  TECHNOLOGY 

A.   INTRODUCTION 

Since,  with  the  exception  of  the  comprex,  relatively  few 
wave  machines  have  been  built,  the  technology  is  not  de- 
veloped to  a  point  where  "preliminary  design"  and  "advanced 
design"  procedures  have  been  formalized.   On  the  contrary, 
each  documented  example  of  a  wave  engine  has  involved  an 
early  arbitrary  decision  on  the  selection  of  rotor  passage 
size  and  orientation  and  an  iterative  construction  of  a  wave 
cycle  that  results  in  a  physically  realizable  system  of 
ports. 

Construction  of  wave  diagrams  for  wave  machines  is  quite 
involved,  requiring  considerable  time  and  effort  for  each 
one.   Complicated  unsteady  wave  phenomena  appear  and  even 
the  simplest  mode  of  operation  requires  calculation  of  an 
array  of  wave  processes.   Approaches  adopted  by  various  in- 
vestigators include  the  method  of  characteristics  and 
graphical  techniques.   Such  "classical"  methods  are  valuable 
in  obtaining  a  proper  understanding  of  the  technology  but 
are,  at  the  same  time,  limited  in  their  ability  to  arrive 
quickly  at  a  preliminary  design  configuration. 

If,  in  the  future,  rational  procedures  are  to  be  devel- 
oped for  selecting  and  optimizing  designs,  it  is  necessary 
to  have  one  set  of  simplified  procedures  to  select  the  most 

16 


desirable  overall  configuration  and  a  more  detailed  analysis 
procedure  with  which  to  finally  optimize.   The  simplified 
procedure  would  be  used  to  select  the  wave  cycle  (and  it  is 
noted  that  passage  stagger  angle  can  be  varied  quite  inde- 
pendently for  a  given  wave  system  to  give  an  infinity  of 
different  machines) ,  and  the  more  detailed  analysis  would  be 
required  to  optimize  the  port  and  passage  designs. 

In  the  present  work,  in  order  to  determine  what  was 
available  and  what  was  needed,  the  task  was  adopted  of  de- 
signing a  wave  rotor  to  perform  as  a  cold-cold  air  "pressure 
equalizer"  (Ref .  6)  .   The  constraints  were  to  use  an  existing 
rotor  with  a  passage  stagger  angle  of  60",  to  use  an  existing 
air  source  at  3  atmospheres  of  pressure,  and  to  pump  air  from 
atmospheric  to  a  common  exhaust  pressure  of  1.5  atmospheres. 
As  a  result  of  the  exercise,  which  will  not  be  documented 
as  such,  it  is  possible  in  the  following  paragraphs  to  pro- 
vide a  review  of  methods  available  for  the  design  and 
analysis  of  wave  rotors.   The  review  includes  a  new  technique 
which  was  implemented  for  preliminary  design  calculations 
(Section  III.B),  progress  made  toward  developing  more  de- 
tailed analytical  approaches  (Section  III.C),  and  the  defi- 
nition of  a  particular  experiment  with  which  to  evaluate, 
in  a  controlled  way,  the  accuracy  of  the  analytical  flow 
descriptions.   The  preliminary  design  technique  and  the 
first  proposed  "turbine  mode"  experiment  are  described  in  a 
recently  submitted  paper  given  in  Appendix  C . 
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B.  PRELIMINARY  CYCLE  CALCULATION  TECHNIQUES 

1.  Riemann  Problem  Solver 

Details  of  a  Riemann  problem  solver  code  developed  at 
NPS  as  a  fast  and  efficient  method  to  carry  out  preliminary 
wave  calculations  are  given  in  Appendix  C,  together  with  an 
example  of  its  application  to  the  proposed  "turbine  mode" 
experiment  wave  cycle.   The  calculational  procedure  requires 
only  minutes  to  carry  through,  with  minimal  requirements  for 
computer  time  and  storage;  e.g.  a  typical  "Riemann  Step"  cal- 
culation requires  1.12  seconds  CPU  time  on  an  IBM  370-3033AP 
computer . 

2.  State  Diagram  Methods 

Unsteady  flows  in  shock  tubes  and  pipes  have  long 
been  treated  by  "classical"  methods  wherein  the  designer  has 
recourse  to  the  construction  of  "state"  diagrams  along  with 
the  wave  diagram.   The  state  diagram  is  a  particularly  useful 
tool  for  wave  rotors  in  that  it  provides  the  designer  with 
some  "feel"  for  the  working  of  the  device.   Based  on  a  refer- 
ence state,  the  state  diagram  is  usually  plotted  in  non- 
dimensional  form,  with  the  ordinate  representing 

Y-1 

2y 
a/a  r-   =  (p/o   -)     and  the  abscissa  representing  u/a   ^. 

ret    '^'^  "  ret  '^  ^  ref 

For  the  case  of  one-dimensional  unsteady  flow  of  an  ideal 
gas,  there  are  two  families  of  physical  characteristics,  the 

C.  and  C   characteristics.   Corresponding  to  the  physical 
characteristics  there  are  two  families  of  state  characteris- 
tics which  relate  the  changes  of  the  flow  speed  and  of  the 
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local  speed  of  sound  along  the  physical  characteristics.   For 

example,  along  the  C_^   physical  characteristic,  whose  path  in 

dx 
the  physical  plane  is  given  by  -^-r-  =  u+  a,  the  values  of  the 

flow  speed  and  the  local  speed  of  sound  Cor  local  pressure 

level)  are  related  by  the  equation  of  the  corresponding  state 

characteristic,  the  T      characteristic,  which  is 

2 
u  +  r  a  =  P  =  constant. 

Y-1 

By  choosing  scales  such  that  a  unit  on  the  ordinate  of  the 

2 
state  plane  is  — y   times  a  unit  on  the  abscissa,  the  state 

characteristics  are  oriented  at  -A5°    in  the  state  plane. 
Following  a  particular  gas  in  a  wave  cycle  from  state  to 
state  is  thus  a  matter  of  following  the  state  characteristics 
in  the  state  plane,  along  the  ±45°  lines.   The  process  can  be 
illustrated  for  the  "turbine  mode"  experiment  described  in 
Appendix  C . 

The  turbine  mode  wave  configuration  is  shown  in  Fig. 
1,  which  is  taken  from  Appendix  C.   The  cycle  is  started  at 
the  origin  of  the  state  plane,  shown  in  Fig.  2,  corresponding 
to  state  1  in  Fig.  1.   The  quiescent  gas  is  accelerated  and 
compressed  to  state  3  by  the  incoming  high  pressure  air. 
This  change  in  state  is  traced  out  simply  on  the  state  plane 
by  moving  along  the  +45"  line  from  the  origin  to  point  3, 
which  is  the  point  of  intersection  of  rhe  state  character- 
istic and  the  incoming  air  "boundary  curve",  drawn  for  the 
specified  inlet  total  pressure  (absolute)  and  wave  rotor 
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geometry.   The  wall  boundary  condition  brings  the  flow  to  a 
stop  with  a  corresponding  increase  in  pressure  (through  the 
reflected  compression  wave) ,  to  state  5  in  Fig.  1.   This 
change  in  state  is  seen  on  the  state  plane  as  line  3-5. 
Assuming,  for  simplicity,  that  the  secondary  waves  in  the 
wave  diagram  (c-e,  e-f,  f-g)  have  a  negligible  influence  on 
the  gas,  states  4  and  5  then  correspond  to  states  7  and  8. 
By  letting  this  high  pressure  gas  out  to  the  outlet  boundary 
conditions  changes  the  state  from  8  to  10  for  the  gas  that  is 
being  tracked.   This  change  is  represented  in  the  state  plane 
as  line  5-10,  which  is  the  state  characteristic  for  the  expan- 
sion wave.   For  the  cycle  to  "close",  the  reflected  expansion 
wave  must  bring  the  gas  back  to  its  original  state.   It  is 
clear  from  the  state  plane  that  this  will  be  the  case  only  if 
the  outlet  conditions  are  at  the  state  defined  by  point  10, 
which  is  the  point  of  intersection  of  line  5-10  and  a  line 
from  the  origin  at  -45°  representing  the  state  characteristic 
for  the  reflected  expansion  wave.   Thus,  for  example,  for 
some  other  inlet  air  conditions,  the  cycle  established  could 
be  as  shown  by  ( 1) - (X) - (Y) - (Z ) - (1)  in  Fig.  2,  with  the  stipu- 
lation that  the  outlet  conditions  are  those  of  point  Z.   Each 
of  such  "closed"  cycles,  however,  give  different  outlet  con- 
ditions, and  due  to  the  strong  interdependence  of  the  parame- 
ters, the  wheel  speed,  output  (shaft)  horsepower,  and  mass 
flow  rate  are  different  for  each  cycle  (Fig.  3). 
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It  is  noted  here  that  the  state  plane  construction 
described  above  is  good  for  small  pressure  ratios  when  only 
weak  shocks  are  involved.   If  two  Mach  lines  of  a  continuous 
compression  wave  intersect  in  the  physical  plane,  then,  as 
long  as  the  pressure  ratio  across  the  resulting  shock  is 
small,  the  straight  state  characteristics  can  be  used.   If, 
however,  the  shock  involved  is  strong,  the  so-called  "shock 
polar"  must  be  used  in  place  of  the  straight  line  characteristic 

Construction  of  state  diagrams  can  become  quite  elabor- 
ate if  supercharger  or  engine  applications  are  involved.   These 
typically  involve  two  different  gases  (combustion  products 
and  fresh  air)  with  differently  oriented  families  of  state 
characteristics  for  each  gas,  if  both  gases  are  tracked  on 
one  state  plane.   Nevertheless,  the  "visual  feel"  that  can  be 
obtained  from  these  diagrams  for  particular  applications  is 
difficult,  if  not  impossible,  to  get  from  purely  numerical 
results.   The  present  day  Brown-Boveri  "comprex"  has  evolved 
over  two  decades  of  research  to  a  high  level  of  sophistica- 
tion as  a  diesel  engine  supercharger.   Even  though  modern 
numerical  codes  are  used  to  calculate  wave  propagation  in 
the  machine,  state  diagrams  are  used  to  assess  the  computed 
results  qualitatively  and  to  monitor  the  changes  which  occur 
from  design  modifications  such  as  the  incorporation  of  pockets 
to  provide  speed/load  range  capability.   The  usefulness  of 
the  state  diagram  approach  coupled  with  modern  numerical  codes 
is  very  evident  in  Ref.  15. 
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3.   Effect  of  Stagger 

Depending  on  the  particular  application  for  which  a 
wave  machine  is  intended,  the  rotor  may  have  axial  passages 
or  "staggered"  passages,  i.e.  canted  at  some  angle  with  re- 
spect to  the  axis  of  the  rotor.   Depending  again  on  the  ap- 
plication of  the  device,  "uniflow"  or  "counterf low"  scavenging 
is  used  corresponding  to  the  gas  (or  gases)  entering  and 
leaving  the  rotor  from  different  ends  or  from  the  same  end. 
The  "turbine  mode"  experiment  described  in  Appendix  C  thus 
has  counter flow  scavenging. 

Uniflow  scavenging  enables  rotor  temperatures  to  be 
minimized  by  bringing  all  parts  of  the  rotor  in  contact  with 
cold  and  hot  fluid  alternately.   This  is  advantageous  for 
engine  applications.   Supercharging  applications  do  not 
usually  encounter  such  high  temperatures  and  may  incorporate 
counterf low  scavenging.   If  the  rotor  passages  are  staggered, 
counterflow  scavenging  makes  work  extraction  simpler  due  to 
the  180°  reversal  that  the  gas  experiences  (thereby  bringing 
about  a  large  change  in  its  absolute  tangential  velocity 
component) .   It  can  be  seen  then  that  the  larger  the  stagger 
angle,  the  greater  is  the  potential  for  work  extraction  for 
the  same  wave  configuration.   The  limiting  factor  here  is 
the  correspondingly  more  complex  port  and  manifold  geometries 
required  for  highly  staggered  rotor  passages.   The  rotor  to 
be  used  in  the  "turbine  mode"  experiment  has  flow  passages 
set  at  60  "^  to  the  axis.   For  the  cycle  described  in  Appendix 
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C,  with  a  mass  flow  rate  m  =  0.754  kg/sec,  rotor  speed 
'^rotor  ~  ^46.2  m/sec,  and  change  in  tangential  velocity 
AV^  =  237  m/sec,  the  ideal  power  that  can  be  extracted  is 
computed  as 

VJ  =  m  •  U   ^    •  AV^ 
rotor    t 

.754  X  146.2  X  237    ^o  ^  u 
=  yyg -  33.6  h.p. 

C.   ENGINEERING  ASPECTS — ANALYTICAL 
1.   Wave  Propagation  Process 

The  wave  propagation  process  is  described  in  some 
detail  in  Appendix  C.   The  emphasis  there  is  on  the  useful- 
ness of  the  Riemann  Problem  Solver  code  for  preliminary  de- 
sign of  wave  rotor  cycles.   Once  such  a  cycle  has  been 
established,  however,  there  is  a  need  for  a  "generalized" 
flow  solver  to  track  the  wave  propagation  and  the  flow  point 
by  point  through  the  rotor  passages.   The  flow  solver  should 
also  incorporate  real  effects  such  as  friction,  heat  trans- 
fer and  leakage  losses . 

Considered  as  unsteady  and  one-dimensional,  the 
governing  equations  for  the  flow  in  the  passages  are  hyper- 
bolic and  lend  themselves  readily  to  solution  by  the  method 
of  characteristics  employing  finite  difference  techniques. 
However,  inherent  in  the  wave  processes  in  real  devices  are 
discontinuities  such  as  shocks  and  contact  surfaces.   The 
detection  of  the  initiation  of  imbedded  shocks  using  the 
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method  of  characteristics  is  based  on  locating  the  point 
where  characteristics  of  the  same  family  intersect.   If  the 
algorithm  used  is  based  on  following  Mach  lines  of  one 
family  and  the  pathlines,  detection  of  initiation  of  shocks 
of  the  other  family  is  difficult.   An  inverse  marching  method 
makes  detection  of  shocks  of  either  family  difficult.   Using 
an  "overall"  direct  marching  algorithm  based  on  following  Mach 
lines  leads  to  fairly  easy  detection  of  initiation  of  shocks, 
but  even  in  this  case,  the  usefulness  of  the  solution  may  be 
sharply  curtailed  if  strong  shocks  are  present  in  the  flow. 
This  is  due  to  the  severe  reduction  in  solution  points  neces- 
sitated by  the  coalescence  of  compression  waves  to  form  strong 
shocks,  leading  to  an  increase  in  the  spacing  between  solution 
points.   The  choice  of  the  grid  and  the  algorithm  used  thus 
depends  strongly  on  the  specific  problem  and  whether  informa- 
tion is  available  beforehand  to  apply  shock-fitting  techniques. 
The  approach  adopted  in  the  present  work  for  numerical 
modelling  of  wave  rotor  applications  uses  the  method  of 
Godunov  to  solve  the  one-  and  two-dimensional  time  dependent 
equations  of  gas  dynamics  in  Eulerian  form.   The  formulation 
is  elaborate,  but  the  method  is  proven  to  treat  any  discon- 
tinuities in  the  flow  (including  contact  surfaces)  more 
realistically  than,  for  example,  artificial  viscosity  methods 
(Ref .  16) .   The  two-dimensional  formulation  is  used  for  those 
regions  of  the  flow  where  such  effects  are  predominant;  for 
example,  rotor  cell  opening  and  closing.   These  results  can 
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then  be  incorporated  to  modify  boundary  conditions  for  the 
overall  one-dimensional  scheme  implemented  to  calculate  the 
entire  cycle.   The  following  section  discusses  some  preliminary 
results  from  a  2-D  modelling  of  the  transient  flow  during  cell 
opening. 

Appendix  D  gives  results  for  the  partial  structure  of 
the  wave  propagation  process  in  the  simple  4-port  cycle  (Ref. 
17)  shown  in  Fig.  Dl.   The  flow  parameters  were  computed  using 
a  new  method  and  new  code  being  developed  currently  at  McDonald- 
Douglas  (Ref.  18),  which  has  previously  demonstrated  an  unusual 
capability  of  following  waves  and  interfaces  from  points  of 
inception,  through  multiple  interactions,  without  distortion. 
As  a  trial,  part  I  (the  high  pressure  part)  of  the  cycle  was 
computed  for  a  pressure  ratio  of  2.0.   Times  t,  and  t2  in 
Appendix  3  refer  to  the  closure  of  the  high  pressure  gas  in- 
let port  and  the  closure  of  the  high  pressure  air  outlet  port 
respectively.   The  shock  wave  generated  by  admission  of  the 
high  pressure  gas  at  time  step  zero  can  easily  be  tracked  for 
successive  time  steps  through  the  pressure  and  velocity  dia- 
grams, while  the  motion  of  the  interface  between  the  com- 
pressed air  and  the  hot  gas  can  be  identified  through  the 
entropy  plots.   The  apparent  "smearing"  of  the  interface 
evident  in  the  plots  is  the  result  of  using  first  order 
(spatial)  accuracy  for  simplicity  during  the  development  of 
the  code,  and  can  easily  be  removed. 
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2 .   Port  Opening  and  Closing 

It  is  suspected  that  strong  2-dimensional  effects  are 
a  major  cause  for  losses  observed  in  the  performance  of  wave 
machines.   These  effects  are  a  consequence  of  finite  cell 
widths  and  the  finite  times  taken  by  the  cells  to  fully  open 
or  fully  close  to  the  various  "stator"  ports.   Although  there 
is  as  yet  insufficient  evidence  in  the  form  of  careful  ex- 
perimental measurements  that  two-dimensional  effects  are  in- 
deed a  major  cause  for  observed  inefficiencies,  it  is  likely 
that  significant  losses  should  be  attributed  to  them. 

The  velocity  and  pressure  contour  plots  presented  in 
Appendix  E  are  preliminary  results  of  tests  toward  a  numerical 
2-D  modelling  of  the  cell  opening  and  closing  process.   The 
code  used  is  an  original  approach  being  developed  in-house  by 
S.  Eidelman*  and  uses  the  2-D  Godunov  method  mentioned  ear- 
lier.  Initial  conditions  for  the  test  case  were  taken  to  be: 

p^  =  5  atm  ,    p^  =  6  kg/m   ,      u^  =  v^  =  0 


)j^  =  1  atm  ,    Pp,  =  1.2  kg/m""  ,    u^ 


The  two  gases  are  originally  on  either  side  of  a  "solid"  wall, 
contained  between  horizontal  parallel  boundaries.   At  the 
start  of  the  integration,  a  diaphragm  is  removed  from  the  gap 
above  the  wall  and  the  wall  is  moved  with  a  constant  velocity 


*NRC  Research  Associate  under  the  NPS  Foundation  Research 
Program. 
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of  200  m/sec  (which  is  representative  of  pitch  line  speeds  of 
wave  rotors) .   The  contour  plots  shown  are  for  seven  different 
time  steps  (corresponding  to  different  position  of  the  wall 
or,  equivalently ,  fraction  of  the  cell  opened),  starting  with 
the  cell  being  open  one-third  of  a  cell  width.   The  time 
periods  shown  are  real  time. 
3.   Ports  and  Pockets 

A  challenging  aspect  of  wave  machine  design  is  in- 
corporating load  and  speed  range  capability  for  off-design 
operation.   Closely  tied  with  this  is  the  question  of  whether 
such  an  operational  range  can  be  obtained  with  a  basic  four- 
port  configuration  typical  of  supercharger  or  engine  appli- 
cations.  No  rigorous  proof  is  offered,  but  it  may  be  shown 
that  serious  practical  difficulties  would  inhibit  the  effi- 
cient working  of  a  four-port  machine  which  has  fixed  port 
geometry,  intermediate  heat  addition  and  from  which  useful 
work  extraction  is  expected.   Appendix  F  outlines  some  of  the 
difficulties,  taking  the  early  Com.prex  design  as  the  test 
configuration. 

Static  "pockets"  in  the  stator  walls  are  possibly 
the  best  design  tool  to  use  for  operational  flexibility  of 
wave  machines.   The  present-day  Comprex  incorporates  three 
types  of  pockets  ("compression",  "expansion"  and  "gas" 
pockets)  to  achieve  speed  and  load  range  capability  (Ref.  15). 
Pearson  used  recirculatory  "loops"  in  his  wave  engine  to 
avoid  losses  inherent  in  pocket  design,  in  effect  adding 
additional  connected  ports. 
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It  is  noted  again  that  state  diagrams  are  an  invalua- 
ble tool  in  visualizing  how  these  elements  help  in  the  opera- 
tion of  wave  machines.   Effort  is  being  made  in  the  present 
work  to  understand  how  best  to  design  and  utilize  "pockets" 
and/or  recirculatory  loops,  and  how  best  to  evaluate  this 
understanding  experixmentally. 

D.   ENGINEERING  ASPECTS — EXPERIMENTAL 

1.  Approach 

A  first  objective  in  the  present  program  was  to  de- 
sign and  conduct  an  experiment  using  existing  hardware 
which  would  serve  to  first  assess  and  subsequently  to  extend 
the  present  level  of  understanding.   The  recommendation  fol- 
lowing the  present  review  is  to  set  up  staged  experiments, 
wherein  unit  processes  of  a  typical  wave  machine,  namely  the 
turbine  cycle  and  the  compressor  cycle,  are  first  modelled 
separately.   Depending  on  the  degree  of  success  achieved  in 
these  tests,  a  combined  cycle  experiment  (which  can  not  be 
simply  a  combination  of  the  two  previous  cycles)  will  follow, 
In  the  longer  term,  "hot"  tests  related  to  a  then  identified 
potential  engine  application  can  be  considered.   The  first 
and  simplest  of  the  proposed  experiments  involves  running 
the  wave  rotor  as  an  impulse  turbine,  as  is  described  in 
Appendix  C, 

2.  Apparatus 

A  schematic  of  the  experiment  is  shown  in  Fig.  4.   A 
detailed  engineering  drawing  of  the  rotor  assembly  is  given 
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in  Appendix  G.   The  drawing  shows  a  half-sectional  elevation 
of  the  assembled  wave  rotor  apparatus  with  all  major  compo- 
nents identified.   Other  partial  views  help  in  visualizing 
the  orientation  of  the  passages  in  the  rotor,  of  which  a 
photograph  is  shown  in  Fig.  5.   The  rotor  and  assembly  hard- 
ware, originally  used  by  the  General  Electric  Company  (Ref. 
4)  is  in  excellent  condition.   Not  shown  in  the  drawing  are 
the  inlet  and  outlet  manifolds  which  must  be  modified  from 
existing  parts.   The  main  constraints  in  the  design  of  the 
inlet  and  outlet  manifolds  are  that  they  must  supply  the 
cycle-dependent  flow  angles  and  peripheral  openings  at  inlet 
and  outlet  ports  and  yet  must  not  interfere  with  other  adja- 
cent components.   (Additional  constraints  would  be  imposed 
in  the  case  of  flight  applications,  such  as  the  low  frontal 
area  required  for  cruise  missile  engines.) 

It  is  noted  that  the  bearing  arrangement  of  one 
cylindrical  roller  bearing  and  a  set  of  duplex  angular  con- 
tact ball  bearings  (assembled  back-to-back)  is  typical  of 
applications  where  it  is  necessary  to  restrict  misalignment 
or  shaft  deflection,  and  where  close  control  of  axial  play 
is  desired.   Clearances  between  the  staror  or  valve  plate 
faces  and  rotor  faces  are  kept  very  small  to  minimize  leakage 
losses.   Any  misalignment  or  axial  movement  could  result  in 
rubbing  of  the  staggered  labyrinths  machined  on  the  rotor  and 
valve  plate  faces. 
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Axial  clearances  can  be  controlled  by  shimming  on 
one  side  of  the  rotor.   This  avoids  having  to  disassemble  the 
entire  machine  in  order  to  set  a  new  clearance.   It  is  ex- 
pected that  the  effect  of  end  clearance,  or  leakage,  will  be 
measured  in  the  course  of  the  study. 

Both  oil  and  forced  air  lubrication/cooling  systems 
may  be  used. 

3 .   Measurement  Program 

As  shown  in  Fig.  4 ,  the  rotor  will  be  attached  to  an 
air  dynamometer  in  order  to  provide^  control  over  the  load 
output  on  the  shaft.   The  load  can  be  varied  over  a  range  of 
the  rotor  speed  by  throttling  the  dynamometer  and  controlling 
the  pressure  ratio  across  the  wave  rotor.   The  "off-design" 
behavior  of  the  rotor  can  be  examined  in  this  way,  and  the 
available  range  and  efficiency  of  operation  established. 

At  the  same  time  the  wave  and  flow  structure  (and 
therefore  wave  timing)  within  the  m.achine  v/ill  be  measured 
using  an  instrumented  closure  plate  on  the  end  remote  from 
the  inlet  and  outlet  manifolds.   By  rotating  the  plate  about 
the  axis,  probes  (pneumatic,  fine-wire  thermocouple  and 
Kulite  semiconductor)  can  be  used  to  peripherally  "survey" 
the  closed  end-wall  in  the  wave  cycle  (Fig.  1) .   A  comparison 
of  the  predicted  and  measured  wave  structures,  pressure  and 
temperature  states  can  be  made. 

Pressures  and  temperatures  are  to  be  measured  at 
locations  within  the  ports  and  in  the  annular  SDace  between 
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the  outer  casing  and  rotor  shroud.   Specific  attention  will 
be  given  to  the  problem  of  measuring  leakage  rates  during 
operation. 

All  data  will  be  acquired  and  reduced  using  the 
Laboratory's  HP  1000  computer  controlled  high  speed  data 
acquisition  system.   It  should  be  possible  eventually  to 
display  on-line  comparisons  of  computed  and  measured  flow 
structure  using  computer  graphics.   It  is  entirely  possible 
that  computer  graphics  could  provide  the  key  to  obtaining  a 
physical  understanding  of  the  on-  and  off-design  behavior 
of  wave  rotors  and  of  the  relative  importance  of  physical 
and  gas  dynamic  changes. 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

To  this  point  in  time,  the  wave  rotor  has  been  an  inter- 
esting concept  which  has  seen  relatively  few  commercially 
successful  applications.   In  most  potential  applications 
however,  a  difficult  development  program  would  have  been 
faced  to  overcome  the  competition  of  an  existing,  or  de- 
veloping simpler  steady-state  alternate  device.   Only  in  the 
case  of  the  Comprex  has  a  competitive  wave-rotor  alternate 
to  a  conventional  turbomachinery  unit  been  successful.   In 
that  particular  case,  as  quoted  by  Rose  (Ref .  7) ,  the  Com- 
prex has  shown  a  performance  which  was  superior  to  a  competi- 
tive standard  turbocharger . 

In  the  case  of  the  Comprex,  at  the  outset,  there  was  a 
potential  for  overcoming  the  slow  acceleration  response 
which  was  characteristic  of  turbochargers .   Without  this 
motivation  the  present  Comprex  would  surely  not  have  been 
developed. 

A  much  stronger  motivation  now  faces  DOD  in  the  area  of 
cruise  missiles.   The  performance  of  small  gas  turbine 
engines  falls  far  short  of  larger  units;  and  how  can  you 
increase  the  cycle  temperature  and  not  have  even  worse  per- 
formance because  of  cooling?   While  the  wave  rotor  may  not 
be  the  perfect  answer,  it  offers  such  a  definite  promise  of 
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advantage  that  its  development  in  an  engine  application 
should  be  seriously  considered. 

The  present  state  of  knowledge  of  wave  rotors  should  be 
seen  realistically.   With  the  exception  of  the  Comprex  there 
have  been  relatively  few  such  devices,  ever.   Individuals 
have  pursued  different  concepts,  largely  in  order  to  simply 
get  them  to  work.   There  is  very  limited  documentation 
available  of  the  wave  rotor  efforts  which  have  been  made,  and 
almost  no  records  of  systematically  obtained  experimental 
data.   On  the  other  hand,  the  work  which  has  been  done, 
notably  by  Pearson,  has  irrefutably  established  the  feasi- 
bility of  developing  a  competitive  small  wave  engine. 

The  most  extensive  bank  of  knowledge  on  wave  rotor  tech- 
nology is  at  Brown-Boveri  Co.  in  Switzerland.   From  the 
paucity  of  design  information  published  on  the  Comprex  after 
twenty  years  of  development  activity,  one  can  conclude  that 
their  design  tools  will  remain  proprietary  to  the  company. 
While  it  may  be  possible  for  DCD  to  contract  a  design  ac- 
tivity to  Brown-Boveri  in  association  with  a  U.S.  engine 
company,  it  can  not  depend  exclusively  on  this  sole  source. 
Furthermore  the  engine  application  has  very  different  con- 
straints and  parameter  requirements  and  it  would  be  unwise 
to  be  limited  to  design  tools  developed  specifically  for 
engine  turbochargers . 

In  this  situation,  overall  it  is  suggested  that  two 
closely  related  programs  are  needed: 


33 


a)  a  prototype  wave  engine  definition  study  for  the  cruise 
missile  application. 

b)  basic  and  applied  research  to  understand,  explore  and 
advance  the  technology  of  wave  rotor  devices. 

The  prototype  engine  definition  study  would  use  current 
knowledge  to  select  the  most  promising  engine  type  and  gener- 
ate a  proposed  hardware  design.   It  is  probable  that  a  rotor- 
component  development  program  would  follow,  and  a  prototype 
engine  program  if  that  was  successful. 

The  technology  research  effort,  to  develop  understanding 
of  wave  rotor  processes,  to  develop  confidence  in  analytical 
and  computational  tools,  to  identify  key  problems  and  explore 
solutions,  is  strongly  recommended  v;hatever  schedule  is 
adopted  in  the  missile  engine  program.   Both  analytical  and 
coupled  experimental  efforts  are  needed.   It  is  emphasized 
that  documented  experimental  experience  is  lacking  most  of 
all,  and  purely  analytical  approaches  can  not  identify  the 
problems  which  need  to  be  solved. 

The  most  technologically  productive  approach  would  be  to 
initially  associate  the  programs  closely,  to  initially  direct 
the  technology  effort  toward  the  cruise  missile  engine,  and 
to,  in  effect,  coordinate  the  effort  (basic  and  applied)  at 
several  centers  toward  the  earliest  realization  of  the  cruise 
missile  application.   This  would  concentrate  resources  toward 
attaining  a  practical  engine,  but  would  also  serve  to  educate 
the  more  academic  investigators  into  the  earliest  realization 
of  the  most  difficult  general  problems  associated  with  the 
technology. 
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It  is  certain  that  this  area  of  technology  has  considera- 
ble potential,  and  computational  and  experimental  tools  are 
available  today  which  should  make  it  feasible  to  understand 
and  therefore  harness  the  unsteady  processes  on  which  the 
technology  is  based. 
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Figure  1.   Simplified  Wave  Diagram  for 

'Impulse  Turbine  Mode'  Operation 
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Figure  5 .   View  of  the  'i'lave  Rotor 
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APPENDIX  A 
SUPPLEMENT  A  OF  REF.  9 

SUPPLEMENT  A 

The  presented  report  is  supplemented  by  two  further  reports  translated 
by  the  Aeronautical  Engine  Laboratories,  Naval  Air  Experimental  Station, 
Naval  Air  Materiel  Center  Philadelphia: 

(a)  Report  on  the  Development  of  a  New  Heat  Energy  Principle  and 
Gas  Turbine  Processes  with  Pressure  Produced  by  Heating  of  the  Working 
Medium.   AEL  Translation  15. 

(b)  The  Heinkel-Hirth  RR  2  (Tuttlingen)  Gas  Turbine  Engine,  Using 
Thermal  Compression.   (Report  by  Heinkel-Hirth,  Stuttgart-Zuf fenhausen, 
December  1946.)   AEL  Translation  56. 

Reports  AEL  15  and  AEL  56  were  not  available  in  their  English  transla- 
tion until  approximately  A  months  after  completion  of  the  main  portion  of 
F-TR-2186-ND.   However,  Dr.  H.  J.  Pabst  von  Chain  consulted  German  text  of 
these  above  referenced  reports  during  the  preparation  of  the  main  report. 
Several  important  theoretical  problems  vital  to  the  thermal  compression 
engine  were  not  treated  in  these  two  Navy  reports. 

The  original  of  AEL  15  was  given  to  the  Navy  in  October  1945,  at  Stutt- 
gart, Germany.  It  discloses  the  general  course  of  development  in  the  field 
of  thermal  compression  by  H.  Wolff,  along  with  a  brief  theoretical  descrip- 
tion of  the  principle. 

After  receiving  AEL  15,  the  U.S.  Navy  ordered  the  construction  of  some 
test  engines  by  Heinkel-Hirth  at  Stuttgart-Zuf fenhausen.   The  test  results 
are  given  in  the  AEL  56.   The  theoretical  analysis  of  the  processes  which 
are  necessary  for  the  evaluation  of  the  test  results  also  are  included  in 
AEL  56. 

Important  theoretical  and  engineering  problems  are  investigated  more 
completely  in  the  T-2  Report,  F-TR-2186-ND,  than  in  the  two  Navy  reports 
AEL  15  and  AEL  56.   F-TR-2186-ND  is  especially  concerned  with  the  following: 

Unsteady  flow  processes 

Comparison  with  the  pressure  exchanger 

Application  for  aircraft 

Suggestions  for  new  test  engines  and  test  program 

More  detailed  description  of  the  method  of  operation 

F-TR-2186-ND,  in  one  volume,  analyzes  the  problems  cited  above;  estab- 
lishes a  technical  history  to  September  1947,  and  presents  a  utilization 
program  for  the  thermal  comcression  engine. 
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APPENDIX  B.   ENGINE  TESTED  3Y  PAUL  W.  HUBER  AT 
NACA  LANGLEY  RESEARCH  CENTER,  CIRCA  19  5  0 
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NOMENCLATURE 

p  =  pressure,  N/m 

3 
p  =  density,  kg/m 

u  =  velocity,  m/sec 

X  =  space  coordinate  (along  rotor  passage) 

t  =  time  coordinate 
H.P.  =  high  pressure 
L.P.  =  low  pressure 

Subscripts: 

1,  2,  3,  .  .  .  ,  10  refer  to  states  shown  in  Fig.  1 

The  terms  'tubes',  'passages'  and  'cells'  have  been  used 
interchangeably  in  the  text. 
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INTRODUCTION 

Wave  rotors  are  devices  in  which  wave  propagation  is  used 
to  effect  a  transfer  of  energy  between  a  gas  and  a  rotating 
shaft  or  directly  between  one  gas  and  another.   Reviews  of 

such  devices  and  their  general  operating  principles  have  ap- 

(1  2) 
peared.   '     Commercial  applications  of  such  devices  have 

been  developed. 

In  the  case  of  direct  transfer  of  energy,  one  gas  (driver) 
at  high  pressure  is  used  to  compress  a  second  gas  (driven) . 
The  process  is  arranged  to  occur  in  tube-like  passages  ar- 
ranged on  the  periphery  of  a  drum,  or  rotor.   The  compression 
is  achieved  by  means  of  compression  waves  or  shock  waves  and 
the  compressed  gas  is  drawn  off  from  the  end  of  the  tube  in 
which  the  process  takes  place.   The  driver  gas  then  undergoes 
a  series  of  expansions  to  a  lower  pressure  and  is  scavenged 
out  by  freshly  inducted  driven  gas  at  approximately  the  same 
pressure  level.   This  fresh  'charge'  is  then  compressed  by 
the  high  pressure  driver  gas  and  the  cycle  repeats  itself. 

Steady  rotation  of  the  drum  sequences  the  passage  of  the 
ends  of  the  tubes  past  stationary  inlet  ports,  cutlet  ports 
and  endwalls.   This  establishes  unsteady  but  repetitive  flow 
processes  within  the  rotating  tubes  and  essentially  steady 
flew  in  the  inlet  and  outlet  ports.   The  tubes  or  passages 
may  be  oriented  axially  or  at  a  'stagger'  angle,  depending 
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on  the  application  of  the  particular  device.   In  general, 
wave  machines  used  as  pure  pressure  exchangers  for  super- 
charging purposes  have  axially  oriented  passages  (e.g.  the 
'Comprex' ,  Ref.  3)  while  those  with  staggered  passages  may 
be  coupled  with  conventional  turbomachinery  or  to  a  drive 
shaft,  since  shaft  work  extraction  is  possible  with  this 
latter  configuration. 

Complicated  unsteady  wave  phenomena  appear  in  the  working 
of  these  devices  and  design  of  even  the  most  simple  mode  of 
operation  requires  calculation  of  an  array  of  wave  processes 
such  as  shock  waves,  reflected  and  'hammer*  shock  waves  (from 
a  wall  or  moving  interface) ,  rarefaction  waves  (connecting 
two  states  of  a  gas  or  produced  by  tube  end  closure) ,  as  well 
as  interactions  between  incident  and  reflected  waves.   In 
general,  the  two  gases  in  these  devices  have  considerably 
different  enthalpies,  leading  to  the  formation  of  contact 
surfaces  which  also  interact  with  the  various  waves.   The 
flow  in  each  tube  is  usually  analyzed  using  one-dimensional 
or  quasi  one-dimensional  approaches,  implying  that  changes  in 
state  occur  only  along  the  passage.   This  enables  the  wave 
processes  to  be  depicted  on  an  x-t  plane.   Such  wave  diagrams 
are  extremely  useful  in  the  examination  of  possible  wave  cy- 
cles, providing  'visual'  information  of  the  wave  paths, 
proper  placing  of  the  ports,  seeing  interactions  and  whether 
a  particular  cycle  'closes':   the  latter  being  necessary  for 
cyclic  operation.   The  construction  of  wave  diagrams  is 
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usually  quite  involved,  requiring  considerable  time  and  ef- 
fort, since  a  mismatch  in  any  one  parameter  gives  rise  to  a 
new  pressure  wave  which  has  to  be  accounted  for  through  the 
complete  rotational  cycle.   The  method  of  characteristics, 
despite  its  limitations  (only  weak  shocks  are  permitted) ,  has 
been  the  technique  generally  used  for  carrying  out  cycle  com- 
putations.     However,  graphical  methods  of  constructing 

wave  diagrams  are  time  consuming  and  may  require  weeks  of 

(4) 
tedious  calculations.      Clearly,  a  fast  and  'user  friendly' 

computational  procedure  is  required  to  carry  out  preliminary 

design  calculations  for  wave  rotor  devices. 

The  method  described  below  offers  a  unified  approach  to 
the  calculation  of  wave  rotor  cycles  with  no  restriction  on 
the  strength  of  the  waves  involved. 

The  method  has  been  exercised  in  the  design  of  a  wave- 
turbine  experiment  at  the  Naval  Postgraduate  School's  Turbo- 
prcpulsion  Laboratory.   The  work  is  part  of  a  program  to 
examine  wave  rotors  and  their  potential  as  components  in 
flight  propulsion  engines. 


METHOD 

The  method  uses  an  approach  described  in  Ref.  5  to  solve 
the  general  ' Riemann  problem' — that  of  finding  the  flow  which 
results  when  two  gases,  each  at  some  specified  initial  state, 
are  suddenly  brought  into  contact  with  each  other.   The  initial 
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states  are  completely  defined  by  specifying  the  pressure, 
density  and  velocity  of  each  gas.   Depending  on  the  initial 
data,  four  types  of  discontinuity  propagations  are  possible, 
involving  combinations  of  shock  and  and  rarefaction  waves 
originating  from  the  point  of  initial  contact  and  proceeding 
to  the  left  or  to  the  right.   The  solution  defining  the  final 
state  is  analytical  and  does  not  impose  restrictions  on  the 
initial  parameters  of  the  gases.   The  type  of  discontinuity 
which  results  is  defined  by  the  solution  and  is  not  prescribed 
beforehand. 

A  computer  code  that  solves  the  Riemann  problem  for  any 
given  initial  states  of  the  two  gases  has  been  developed. 
(The  program  in  PASCAL  is  available  from  the  authors  upon 
request.)   The  program  is  set  up  such  that  the  gas  with  the 
higher  initial  pressure  may  be  on  either  side.   The  results 
give  the  type  of  discontinuity,  e.g.  shock-shock,  shock- 
rarefaction,  etc.,  the  pressure  and  velocity  at  the  inter- 
face and  the  densities  on  either  side  of  the  interface. 
Velocities  of  propagation  of  the  waves  involved  are  also 
computed,  with  two  velocities  calculated  for  rarefaction  fans 
corresponding  to  the  head  and  tail  waves.   The  pressure  at 
the  interface  is  computed  by  an  iterative  procedure,  using 
the  average  of  the  two  initial  pressures  as  the  initial  value. 
Entropy  checks  are  incorporated  into  the  program  which  dis- 
allow the  appearance  of  expansion  shocks.   Although  the  pro- 
gram, in  its  present  form,  can  handle  only  gases  which  have 
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the  same  ratio  of  specific  heats,  it  can  quite  readily  be 
modified  to  deal  with  different  specific  heat  ratios,  as  is 
outlined  in  Ref.  6.   The  proposed  wave-turbine  and  subsequent 
experiments  involve  only  relatively  cold  air  without  combus- 
tion products,  and  hence  cycle  calculations  with  different 
specific  heat  ratios  have  not  been  necessary. 

The  following  section  describes  the  wave-turbine  as  an 
example  of  how  cycle  calculations  may  be  carried  out  using 
the  program. 


EXAMPLE 

Figure  1  shows  the  wave  diagram  (not  to  scale)  for  the 
proposed  'turbine  mode'  experiment  in  which  a  wave  rotor 
would  be  used  to  work  as  an  impulse  turbine  and  produce  shaft 
power  output.   For  clarity,  the  waves  have  been  shown  as 

single,  straight  lines,  with  full  lines  indicating  compres- 
sion or  shock  waves,  chain-dotted  lines  indicating  contact 
discontinuities  and  dashed  lines  representing  expansion  fans. 
The  direction  of  rotation  of  the  rotor  is  as  indicated  by  the 
arrow  at  the  bottom  of  the  figure.   The  right  side  of  the 
rotor  is  blocked  off,  and  the  left  side  tube  ends  open  to 
inlet  and  outlet  ports  alternately.   Ti^a  encircled  numbers 
depict  regions  of  uniform  flow  but  at  different  state  with 
respect  to  adjoining  regions.   Subscripted  flow  parameters 
used  in  the  following  discussion  correspond  to  these  states. 
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starting  the  cycle  at  state  1  at  the  bottom  of  the  fig- 
ure, the  rotor  tubes  are  filled  with  air  in  a  quiescent  state 
and  approximately  ambient  pressure.   The  cells  arc  then 
brought  into  contact  with  incoming  high  pressure  air  at  the 
inlet  port.   This  generates  a  shock  wave  Ca-b)  which  propa- 
gates into  the  air  at  state  1,  raising  its  pressure  and  den- 
sity to  that  of  state  3.   A  mass  velocity  u^  is  also  generated 
behind  the  shock.   An  interface  (a-c)  separates  the  incoming 
air  and  the  compressed  air,  although  for  this  'cold'  con- 
figuration the  two  densities  are  not  very  different. 

Shock  (a-b)  reflects  off  the  solid  boundary  at  b  as  (c-b) 
and  intersects  the  slower  moving  interface  at  point  c,  where 
part  of  the  shock  is  transmitted  (c-d)  and  part  is  reflected 
(c-e) ,  creating  states  4,  5  and  6.   The  velocity  in  state  5 
is  zero  and  nearly  so  in  states  4  and  6.   The  inlet  port  is 
closed  when  the  transmitted  shock  (c-d)  arrives  at  the  left 
end.   Shock  (c-e)  gets  reflected  again  at  point  e  as  (e-f) 
and  continues  on  as  (f-g) ;  these  'secondary'  reflections  are 
of  almost  zero  strength  and  do  not  affect  the  flow  significantly, 

At  point  g,  the  air  outlet  port  is  opened  and  the  compressed 
air  in  state  7  exits,  with  an  expansion  fan  propagating  in  the 
opposite  direction.   The  interface  bends  towards  the  exit  and 
reaches  the  end  of  the  passage  at  point  j .   The  arrival  of 
the  interface  at  point  j  is  timed  to  match  the  arrival  of  the 
reflected  expansion  wave  (i-j),  and  the  outlet  port  is  closed 
at  this  moment.   The  air  in  the  cells  is  again  in  a  quiescent 
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state  and  should  be  at  the  same  pressure  and  density  as  that 
of  the  original  state  1  for  the  cycle  to  be  repeated. 

The  passages  in  the  rotor  are  'staggered'  at  60°  to  the 
axis  of  the  rotor,  and  the  change  in  the  angular  momentum 
effected  by  the  reversal  in  direction  of  the  air  flow  pro- 
vides a  torque  and  allows  the  extraction  of  shaft  work. 

The  calculations  using  the  Riemann  solver  code  are  as 
follows : 

STEP  1;   H.P.  air  at  state  2  hitting  L.P.  air  at  state  1. 
Initial  conditions  for  Riemann  problem  solver: 


P9/  ^9   Pi  /  P-1  /  ^1  =  0 


2'    ^2'       2 


Discontinuity  Type:   Rarefaction-Shock 
The  shock  wave  propagates  to  the  right  into  the  rotor  passage 

and  the  rarefaction  propagates  to  the  left  into  the  inlet 
port.   A  contact  surface  follows  the  shock  wave  into  the  pas- 
sage at  a  slower  speed.   In  the  limiting  case  where  no  dis- 
turbances are  to  be  propagated  into  the  port  (for  uniform 
flow  conditions) ,  the  initial  parameters  may  be  varied  to 
obtain  the  singular  condition  in  which  no  wave  propagates  to 
the  left  and  a  shock  wave  propagates  to  the  right.   This  can 
be  described  as  the  discontinuity  type  'no  wave-shock'. 

State  3  is  then  completely  defined. 
STEP  2:   The  flow  in  region  3  is  confronted  by  the  'wall' 
boundary  condition  at  the  right  side  and  is  brought  to  a  halt 
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by  the  reflected  shock  (b-c) .   This  situation  is  analogous 
to  the  flow  colliding  with  a  flow  at  the  same  state  but  with 
equal  and  opposite  velocity. 
Initial  Conditions: 


P3'  P3'  "^3      P3'  P3'  "'': 


Discontinuity  Type:   Shock-Shock 

State  5  is  completely  defined. 
STEP  3 ;   The  reflected  shock  (b-c)  intersects  the  interface 
at  point  c. 

Initial  Conditions: 


P2/  P2'  ^2   ^5'  ^5'  ^5  =  ° 


Discontinuity  Type:   Shock-Shock  (or  Shock-Rarefaction 

depending  on  incoming  flow 
temperature) 

States  4  and  6  are  completely  defined. 
STEP  4 :   Reflected  shock  (c-e)  is  re-reflected  from  the  wall 
at  point  e. 

Initial  Conditions: 


P6'  P6'  ""e  ^6'    ^6'  "^i 


Discontinuity  Type:   Shock-Shock 
State  8  is  completely  defined. 
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STEP  5;   Shock  (e-f)  hits  interface  at  point  f. 
Initial  Conditions: 


^6'  ^6'  ^6   P6'  P6'  -^( 


Discontinuity  Type:   Shock-Rarefaction 

State  7  is  completely  defined. 
STEP  6 ;   Air  in  the  rotor  cells  is  released  to  exhaust  (am- 
bient) conditions. 

Initial  Conditions: 


STEP  7 


Pamb'  Pamb'  ^amb  =  ^   ^7'  ^1'    ^7  =  ° 


Discontinuity  Type:   Rarefaction-Shock 

State  9  is  completely  defined. 

Rarefaction  wave  (g-h)  intersects  interface  at  point  h, 

Initial  Conditions: 


Pg/  Pg/  -^9   P9/  Pgr  ^8  =  ^ 


Discontinuity  Type:   No  wave-Rarefaction 

States  9  and  10  are  defined.   State  9  should  be  the 

same  as  the  one  obtained  previously  in  Step  6. 

STEP  3:   Rarefaction  wave  hits  the  wall  at  point  i,  and  is 

reflected  in  like  sense. 
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Initial  Conditions: 


'10'  ''lO'    10   -10'  -"10'  "^10 


Discontinuity  Type:   Rarefaction-Rarefaction 
State  1'  is  defined.   This  should  match  with  the 
original  state  1  for  'cycle  closure',  which  is  required  if 
the  rotor  is  to  operate  continuously. 

Table  I  gives  the  values  computed  for  the  cycle  and  with 
the  procedure  described  above.   The  pressures  are  static  values 
and  the  velocities  are  referred  to  the  rotor. 

Clearly,  a  mismatch  of  rotor  speed  and/or  inlet  conditions 
from  design  point  values  will  cause  new  waves  to  be  generated 
which  have  to  be  carried  through  the  entire  cycle  in  order  to 
assess  their  overall  effect  on  the  performance.   (This  is  done 
in  practice  to  incorporate  features  such  as  'pockets'  which 
make  these  devices  operationally  flexible,  see  Refs.  1  and  3.) 
For  preliminary  design  purposes,  however,  having  to  use  a 
generalized  flow  solver  to  arrive  at  a  preliminary  wave  dia- 
gram would  be  time  consuming  and  expensive.   The  Riemann  pro- 
gram is  particularly  useful  for  this  purpose  because  of  its 
'building  block'  approach,  which  allows  walking  through  any- 
wave  configuration  state  by  state.   Once  a  viable  cycle  is 
established,  a  detailed  flow  solver  can  be  applied  to  incor- 
porate effects  of  friction,  heat  transfer  and  the  finite  times 
taken  for  cell  opening  and  closing. 


56 


fH 
O 

u 


<0 

c 

•H 

Eh 

(U 
Ul 

rH 

3 

S 


0 
vw 

to 

QJ 

3 

.-H 

> 

4J 


s 
o 
u 


CTN 

vO 

vO 

o 

» 

m 

•* 

o 

iH 

>3- 

• 

ON 

rH 

-H 

«N 

m 

t— I 

i-H 

m 

• 

O 

CSI 

1— 1 

r^ 

rH 

• 

• 

<m 

P- 

CN 

rH 

rH 

1 

^ 

00 

00 

^ 

-d- 

iri 

1^ 

• 

a\ 

1— i 

i—t 

SO 

r^ 

• 

cn 

. 

CN 

-H 

.-t 

1 

vO 

m 

r~« 

O 

LO 

00 

00 

I— t 

o 

o 

• 

li 

• 

m 

CM 

1^ 

CO 

rH 

1—1 

CN 

1^ 

p-» 

CN 

o 

ON 

• 

Jl 

• 

f^ 

CN 

<r 

ON 

rH 

^ 

in 

vO 

r-~ 

1-H 

o 

ON 

• 

l\ 

• 

fO 

OJ 

ON 

CO 

m 

ro 

VT 

m 

vO 

<-i 

o 

ON 

• 

en 

CN 

-T 

ON 

rH 

nH 

CN 

<r 

r^ 

CN 

o 

ON 

• 

m 

<VJ 

c^ 

LO 

00 

00 

-T 

m 

• 

n 

r^ 

fH 

vO 

r^ 

• 

n 

• 

CN 

-H 

— t 

m 

m 

m 

oo 

■<r 

00 

• 

CNJ 

rW 

.-! 

vO 

r^ 

• 

ro 

CN 

1—1 

rH 

00 

^ 

00 

ON 

00 

tH 

ro 

<r^ 

o 

■^ 

* 

ON 

-H 

1 

c3       y 

r-s 

o         / 

CN 

r-t            / 

2 

ao      / 

""^^ 

<u     / 

0)          3 

>% 

05     / 

M 

>.          -^ 

U                   ^— V 

3       'J-1 

-J       c-^ 

•H              CJ 

W         1 

•H              5 

•J             CU 

cn        O 

en       -^ 

C            M 

OJ         -^ 

a        00 

f— '           ^^         j 

u        X 

OJ           J^ 

a>         S 

CU            "—' 

Q         ^--' 

>        ^^      1 

57 


CONCLUSIONS 

The  calculational  procedure  in  the  example  above  required 
only  minutes  to  carry  through,  with  minimal  requirements  for 
computer  time  or  storage.   (On  average,  a  typical  'Riemann 
Step'  calculation  required  1.12  seconds  CPU  time  on  an  IBM  370* 
30  3  3AP  computer.)   The  Riemann  problem  solver  code  therefore 
gives  a  fast,  efficient  and  unified  approach  to  carry  out  the 
preliminary  design  of  wave  rotor  devices  with  diverse  wave 
structures  and  pressure  ratios. 

It  is  noted  that  the  code  may  be  coded  easily  on  any 
'home'  computer  and  does  not  require  external  hardware  or 
software  libraries. 


CQ 


L.P.  Air 
Outlet 


H.P.  Air 
Inlet 


Figure  1.   Simplified  Wave  Diagram  for 

' Impulse  Turbine  Mode '  Operation 
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APPENDIX    D 


PRELIMINARY    COMPUTATIONAL    PJISULTS 
FOR    SIJ^PLE    4-oOR'^    WAVE    riACIIINE 
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APPENDIX  E 

PRELIMINARY  COMPUTATIONAL  RESULTS 
FOR  CELL  OPENING/CLOSING  PROCESS 
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APPENDIX  F 
THE  Si:4PLE  FOUR-PORT  DIRECT-FLOW  ROTOR 

The  simplest  arrangement  of  ports,  waves  and  interfaces 

which  can  be  used  in  a  direct  flow  wave  rotor  is  shown  in  Fig. 

Fl.   The  process  is  illustrated  on  a  T-S  diagram  in  Fig.  F2 

and  the  velocity  diagrams  at  the  ports  are  illustrated  in  Fig. 

F3. 

The  states  shown  at  stations  Q^  and  Q   in  Fig.  Fl  are  qui- 

a      g 

escent  relative  to  the  rotor,  the  relative  flow  velocities 
W-,  =  W.  and  W„  =  W^ .   From  continuity  across  Q  ,  the  flow  rate 
of  air,  w  ,  is  given  by 

w   =  D  ,'V  -N  F(l) 

a    ■ qa   c   u 

where  o    is  the  density  of  the  air,  V  is  the  cell  volume  and 
•  qa  ^  '   c 

N   is  the  rotor  velocity  excressed  as  cells  oer  unit  time. 

Similarly,  continuity  across  Q   requires  that  the  flow  rate 

y 

of  gas,  w  ,  is  given  by 
9 

w   =  p   'V  -N  F(2) 

g    ^qg   c   u 

where  p    is  the  density  of  the  quiescent  gas  in  the  passages. 

From  Eq.  F(l)  and  Eq.  F(2), 

w     p 

^  =  ^^  F(3) 

^g   'qg 

Since  the  pressures  at  stations  Q   and  Q   are  of  compara- 

a      g  '^ 

ble  magnitudes,  differences  in  temperature  from  the  heat 


p 
addition  will  cause  the  ratio  —2—  to  become  larger  than  unity. 

■  qg. 

Equation  F(3)  then  requires  that  w  >w   so  that  not  all  the 
^  ^  a    g 

air  which  is  compressed  can  be  passed  back  through  the  rotor. 
The  rate  of  work  which  can  be  obtained  from  the  shaft  can 
be  calculated  from  momentum  considerations.   Since  the  flows 
at  the  ports  are  steady,  the  axial  moment  exerted  on  the  shaft 
in  the  direction  of  rotation,  iM  ,  is  given  by: 

where  R  is  the  radial  distance,  so  that  if  oj  is  the  rotation 
rate  the  rate  of  work  output  on  the  shaft,  W   is  given  by 

^o  "  -\^  "  ^'"^^^1  "  "^U2^'^a  ^  ^"^^'-^3  "  ^'^^4^^g      ^^^^ 
In  general,  from  Fig.  F3, 

V   =  W  sin  6  +  U  F(6) 

u 

so  that 

W   =  (UW,sin  3  -  UW^sin  3)w,  +  (UW^sin  3  -  UW  .sin  B)  w    F(7) 
0I2       a      3         4       g 

Since  W .  =  W,  and  W2=W^, 

W^  =  U  sin  3(W^  -  W^)  (w^  -  w  ) 

w  j 
=  ;V  |1  -  -7^  (W,  -  W^)Usxn3  F(8) 

a     Wa  !   -'-    "^ 
L      ^1 

It  is  interesting  to  note  that  no  work  can  be  obtained  if 

the  rotor  cells  are  axial  (3=0),  but  that  work  can  be  extracted 


*Vavra,  M.  H.,  "Aero-Thermodynamics  and  Flow  in  Turbo- 
machines,"  John  Wiley  &  Sons,  New  York,  19  60. 
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with  "straight"  passages  if  they  are  staggered  (2  j^  0).   Equa- 
tion F(3)  shows  that  as  heat  is  added  (p    becominc  smaller 

qg 

than  p   ) ,  it  is  required  that  -7^  become  progressively  smaller 
ga  ^a 

than  unity.   From  Sq.  FC3),  progressively  more  work  can  then 
be  extracted  on  the  shaft,  with  the  stipulation  that  the  low 
pressure  scavenging  velocity  W,  be  greater  than  the  high  pres- 
sure scavenging  velocity  W^ . 

For  axial  passages,  when  the  nett  work  output  is  zero  from 
Eq.  F(8),  the  energy  from  the  gas  is  transferred  directly  to 
the  air.   The  steady-flow  energy  equation  then  gives 

w,  Cp^(Tt2  -  Tti)  =  Wg  Cpg(Tt3  -  Tt4)  F(9) 


a 


or 


■"t3  -  ■^■t4i    cp^ 


T4-^  -  Ti-J    c. 


^Tt2  -  TtiJ    Cpg   ^wg_ 


^a 


F(10) 


if  the  efficiencies  of  compression  and  expansion  are  de- 
fined in  the  usual  way;  referring  to  Fig.  F2 , 

^t2is  -  ^ti 
^t-t     ^t2    ^^1 

is  the  total-to-total  efficiency  of  compression  and 

Tt3  -  Tt4 

"^T^    =  FS ^ —  F(12) 

Tt-s    Tt3  -  T4.^ 

is  the  total-to-static  efficiency  of  compression. 

Using  Eq,  F(ll)  and  Eq.  F(12),  Sq.  F(10)  can  be  rewritten 
as 
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t3 


1  - 


T4  . 


•  r^'^ 


(T 


t2 


IS 


tl   T 


tl 


-  1 


l=Pg. 


w 


w. 


If  it  is  assumed  that  pt3  =  p^^  ^-^^  that  p, 
total  pressure  ratio  of  the  cycle  is  defined  as 


F(131. 


pt-j,  and  the 


Pt 


n    = 


Pt3 


Pti 


F(14) 


then  Sq.  F(13)  can  be  written  as 


_t3 
'tl 


Y-1 


'^C+.-i-  •  ^Tt-- 


t-s 


'Pa 


'P 


qj 


w 


w. 


F(15) 


Since  the  air  to  gas  weight  flow  ratio  is  controlled  by 
the  heat  added  through  Sq.  F(3)  the  expression  in  Eq.  F(15) 
relates  the  heat  which  can  be  added  in  such  a  cycle  to  the 
pressure  ratio  and  energy  transfer  efficiencies.   However, 
it  is  possible  that  in  practice,  Eq.  F(15)  can  not  be  satis- 
fied for  pressure  ratios  and  efficiencies  which  are  realistic 

The  particular  case  of  a  4-port  rotor  with  axial  passages 
has  been  considered  in  an  approximate  way  by  Hussmann  (Ref. 
19)  and  more  recently  by  Berchtold  (Ref.  17).   Hussmann  has 
also  discussed,  somewhat  empirically,  the  apparent  inability 
of  the  simple  4-port  machine  to  maintain  proper  scavenging 
without  recourse  to  external  aids  such  as  blowers  or  fans. 
It  is  possible  that  for  effective  operation  with  good  sca- 
venging, both  outgoing  velocities  (and  especially,  the  low 
pressure  scavenge  velocity)  may  have  to  be  increased. 
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->  Wa,r  -    W^o^ 


p-      —      ^      —      HEAT     AODl  T»ONi 


_      ^ 
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w, 
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•  w 
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-■> 
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1                       ^ 
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